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EFFECTS OF CLOSURE ON MAJOR AND TRACE ELEMENT VARIATIONS IN APOLLO 11, 12 
AND 15 MARE BASALTS AND REGOLITH SAMPLES 


ABSTRACT 

Major element and major element plus trace element analyses have been 
selected from the lunar data base for Apollo 11, 12 and 15 basalt and regollth 
samples. Summary statistics for each of the six data sets have been compiled 
and the effects of closure on the Pearson product moment correlation coeffi- 
cient have been investigated using the Chayes and Kruskal approximation proce- 
dure. In general there are two types of closure effects evident in these data 
sets: (1) negative correlations of Intermediate size which are solely the 

result of closure and (2) correlations of small absolute value which depart 
significantly from their expected closure correlations which are of intermed- 
iate size. The Apollo 15 regolith samples do not conform to the Chayes and 
Kruskal model because Si02 is behaving as a constant and this set of analyses 
can not be meaningfully compared with the other five data sets. Particular 
attention was directed towards the possibility of inducing a strong positive 
correlation between variables present in trace amounts as a result of closure. 

It is shown that a positive closure correlation will arise only when the pro- 
duct of the coefficients of variation is very small (less than 0.01 for most data 
sets) and, in general, trace elements in the lunar data sets exhibit relatively 
large coefficients of variation. Preliminary experiments suggest that closure 
strongly effects both the Spearman and Kendall rank correlation coefficients 
and, at the present time, the closure effect can not be accounted for. These 
rank correlations may be better suited to assessing the strength of associa- 
tion between geochemical variables than is the Pearson product moment corre- 
lation coefficient and simulation experiments are being continued. 
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EFFECTS OF CLOSUBE ON MAJOR AND TRACE ELEMENT VARIATIONS IN APOLLO 11, 12 
AND 15 MARE BASALTS AND REGOLITH SAMPLES 


Introduction 

It is a well established fact that percentages fanned from uncorrelated 
parent variables will be correlated as a result of percentage formation 
(summarized in Chayes, 1971). As an extreme examples, one can form percen- 
tages from a set of binary data (% X + % Y = 100%), and regardless of the 
correlation between X and Y the correlation between X and Y will be -1.0. 

If the parent variables are homogeneous in mean and variance the correlation 
Induced by percentage formation Is given by: 

p^j “ -1,0/(M-1.0) Equation (1) (Chayes, 1971) 

where M is the number of variables in the closed data set (set of percentages) 
As M increases the magnitude of the Induced correlation decreases (p^^^ ■ -0,5 
for M “ 3 and p^^ ■ -0.1 for M * 11). However, a data array of geochemical 
interest with variables homogeneous in mean and variance is extremely unlikely 
so that one can not argue that the closure effect can be reduced simply by 
Increasing the number of variables in the closed data array. 

Chayes and Kruskal (1966) developed an approximation n«thod for computing 
the correlation that will arise solely as the result of percentage formation 
from uncorrelated parent variables. This closure correlation (P^j) can be 
used as the null value for testing an observed correlation for significant 
departures from randomness. Indlscrimlnant use of zero as a null value or 
a visual interpretation of the strength of a variation diagram may result in 
believing that a correlation is significant when in fact closure alone is 
responsible. Similarly, correlations of relatively small magnitude, which 
require accepting the null hypothesis when zero is used as the null, may be 
significant if the closure correlation is relatively large. 
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Whole sample major and trace element analyses of Apollo 11, 12 and 
15 basalts and regollth samples have been compiled and closure effects 
assessed. Particular attention has been directed towards examining the 
possibilities of induced positive correlation between variables present 
In trace amounts. 

The Chayes and Kruskal Test 

Approximation of the expected closure correlations for the observed 
data array U requires estimation of the means and variances of the variables 
in the X array - the hypothetical closed array - which is characterized by 
possessing zero covariances and hence zero correlations. When X is closed 
to form Y - the hypothetical closed array - the means and variances of the 
variables in Y are equal to those in the observed data array U. As the 
correlation in X are zero, the correlations in Y are solely the result of 
closure. 

The means of the variables in X are the observed mean proportions of 

the variables in U (Chayes, 1971). Hypothetical open variances can be 

computed from the observed mean proportions and variances of the parent 
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variables in U. The sum of the open variances (o^ ) is given by: 


M 


0 ^ • 1-1 


(s^V(1.0-2.0y)) 


M 


1-1 


Equation (2) (Chayes, 1971, 
Eq. 9.1) 


((yj^(1.0-y^))/(1.0-2.0y^)) 


where y^ and s^ are the observed mean proportion and variance of the ith 
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variable respectively. Each open variance (o^ ) can be computed from: 


-2 

2 ^1 


- °t^> 


(1.0-2.0y^) 


Equation (3) (Chayes, 1971, 
Eq. 9.2) 
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where Is the observed coefficient of variation of the 1th variable. 

If none of the variables has a mean proportion greater than O.SO an 

2 2 

open variance will be negative If and only If Is less than 
(Equation 3). 

Assuming that all of the open variances are positive one can compute 

the correlation Induced by closure (pj^j) as follows: 

“ - 2 - 2 - 2 

- y Oj^ - y^Oj )/Sj^8j Equation (4) (Chayes, 1971, 

Eq • 6,1) 

- 2 

where y , s and 6 are the observed mean proportion, the observed standard 
■'n n n r r . 

deviation and the open variance of the nth variable respectively. 

In general, the greater the variances of the variables the greater the 
induced negative cor-elatlon and, for variables with small variances, the 
induced correlation may be positive. Chayes and Kruskal (1966) noted that the 
conditions for positive null correlation would most likely be satisfied by 
variables with small means and variances and concluded that null correlations 
between variables present In trace amounts would often (If not always) be 
positive. 

Thus, the investigator Is faced with a definite problem in assessing 

the strength of a geochemical variation diagram and/or in testing an observed 

correlation for significant departures from randomness. The aprlorl use of 

zero as a null value must be avoided for percentage data and. In general, 

each correlation must be considered to possess its own null value. Also, 
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one must have a complete chemical analysis to evaluate o^ ; that is, partial 
chemical analyses can not be tested using the Chayes and Kruskal model. 

Effect of Closure on Correlation Between Major Elements 

The lunar data base (updated version of March, 1976, supplied by J. 
Warner of JSC) served as the primary data base. Subsets of Apollo 11, 12 
and 15 basalts and regollths have been extracted at the present time. Each 
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subset has been carefully examined, duplicate analysis or entries removed 
and anomalous entries verified by checking the original reference. These 
corrected subsets are serving as the data bases for continuing studies of 
vithin and between mission natural groupings (Butler, 1976), theoretical 
frequency distributions of major and trace elements and the effects of 
closure* 

Vfhole sample (coded ALL in the lunar data base) major element and 
major plus trace element analyses were selected for each mission using the 
following criteria: 

(1) Si02, ^120^, Ti02, MnO, MgO, CaO, Na20, K^O and Cr202 were analyzed 
for by the same investigator (s) and published in the san« reference; 

(2) the sum of the major elements was between 98.5 and 101.5 

(3) trace element analyses were published in the same original reference 
as were the major element analyses. 

Originally it was desired to Include P 2®5 ^ apparent 

paucity of information concerning these two variables (especially in the 
Apollo 15 basalt and regolith samples). Whenever possible S and ^ 2^5 
analyses or averages were included to provide more complete data arrays. 

In general there are approximately 2 to 3 times as many major element 
analyses as there are major plus trace element analyses for a given mission. 
In fact, there were less than 10 samples of majo plus trace element analyses 
for the Apollo 11 and 12 regolith samples so these were excluded from fur- 
ther consideration. 

Several previous investigators have published averages of basalts and 
regollths (for example, Taylor, 1975 and Rose et al., 1976). Frequently 
such averages have been computed using limited portions of the total 
available information and I am not aware of published averages making use 
of all available information. This is especially critical for least squares 
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mixing model studies. However, one must question the meaning of any set 

of computed averages. For example, there ore some samples for which many 

analyses are available (15555 for example) and averages created without 

weighting may prove difficult to interpret. If one computed an average 

composition for each sample and averaged the results to give an average 

composition for the mission the weighting would be equal for each sample 

and not enough Is known about the sampling efficiency at each site to 

warrant an equal weighting assumption. Difficulties arise, for example, 

when averages of basalt and regollth samples are compared. It one has 

confidence In the averages ne may be compelled to postulate the presence 

of an "exotic" component or process to account for any observed differences. 

Taylor (1975) concluded that the "true abundance" of aluminous basalts has 

been underestimated at the Apollo 11 and 12 sites with the result that a 
source for the higher AI 2 O 2 of the regollth must be identified. 

Summary statistics for the major element analyses which satisfied the 

prevloxisly cited criteria are given in Table 1 (mean, coefficient of variation 
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(C) and % of the total variance (% s )). More detailed sets of summary 
statistics for each analyzed sample may be obtained by writing the author. 

Use of such summaries, however, must be undertaken with considerable discretion. 

Chayes (1965) advocated computing the sum of the variances as a measure 
of the total variability of the data array. From Table 1 it is evident that 
the regollth samples have significantly lower total variances than tha basalts 
for a given mission. Also, the variability of a given variable (as estimated 
from the coefficient of variation) appears to be greater for the basalts than 
for the regollths for a given mission. In part this might be the result of 
fewer samples in the regollth suite but most likely reflects the "gardening" 
processes accompanying regollth formation. 
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Open vatlances were computed for each of the 6 data sets sutranarized 

in Table 1. The open variance for SIO 2 In the Apollo 15 regollth suite 

Is negative so that individual correlations can not be tested for significant 

departures from randomness. The coefficient of variation of SiO^ is 0.013 
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for the Apollo 15 regollth samples and C Si02 is less than the sum of the 

open variances (0.00456) so that the open variance for Si02 is negative 
because SIO 2 is behaving as a quasi constant. All of the remaining 5 suites 
possess positive open variances. Rather than list all 66 correlations for 
each of the 5 data sets that conform to the Ciiayes and Kruskal model it was 
decided to list only those for which the result of testing against the clo- 
oure correlation (Equation 4) significantly differed from testing against 
a null of zero (Table 2). 

It Is evident from Table 2 that two types ot misinformation (butler, 
1976) are present in the 5 sets of lunar chemical analyses: 

(1) correlations which would be significant If tested against a null 
of zero but fail to significantly depart from their expected 
closure correlation; 

(2) correlations of small absolute value which would require acceptance 
of the null hypothesis If tested against zero but significantly 
depart from their intermediate in size closure correlation. 

It is apparent from Table 2 that large negative closure correlations occur 
between variables with large variance. Of particular interest are the 
relatively large positive closure correlations which will be considered in 
the following section. 

Many geochemists have been using multivariate statistical methods which 
require that the matrix of correlation coefficients be used as a measure of 
similarity for an R mode analysis. As demonstrated previously (Butler, 1976) 
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such methods can produce biased results unless the effects of closure are 
somehow taken Into account. In Tables 3, 4« and 5 correlations are given 
for Apollo 11, 12 and 15 respectively. Those that are significant at the 992 
confidence level when tested against their appropriate closure correlations 
are Indicated by either a -t- or - sign. Correlations for basalts are given 
in the upper half and correlations for regoliths are given in the lower 
half of the matricies in Tables 3, 4, and 5. There Is a great deal of infor- 
mation contained in Tables 3, 4 and 5 and more detailed comparlslons are 
currently being prepared. However, it is apparent that there are major 
correlation coefficient pattern differences between basalt and regolith 
sacq)les from the Apollo 11 and 12 sites. It Is likely (as many previous 
Investigators have noted) that the regolith at a given site Is not derived 
by simple mechanical breakup of basalt (s) with an average composl-lon of 
that given In Table 1. 

As noted previously, the Apollo 15 regolith samples do not conform to the 
Cltayes and Kruskal model. Suggested stratlgles for the elimination of nega- 
tive open variances (Chayes, 1971 and Butler, 1975) have involved using 
normative quantities of the alkalis (KAIO 2 , orthoclase or albite). Such 
transformations, however, were designed for the elimination of negative 
open variances for AI 2 O 2 and Ns 20 in data sets in which SIO 2 was greater 
than 50%. Rational transformations for removing the negative open variance 
of Si02 in the Apollo 15 regolith suite and in the three postulated suites 
of Apollo 15 basalts (Butler, 1976) are currently being investigated. 

The Effect of Closure on Correlations Between Trace Elements 

Trace elements possess very small variances as compared with the variances 
of the major elements and one would not expect that the closure correlations 
between trace elements would be negative and of Intermediate size. Positive 


correlation Is possible If and only If the following (a rearrangement 
of Equation 4) holds: 
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Equation (5). 


As a complete chemical analysis Is required to compvjte , trace element data 
sets were compiled from the sets of major element analyses summarized In 
Table 1. Relatively few major plus trace element analyses have been pub- 
lished so an attempt was made to maximize the number of samples In each of 
the four. selected arrays (Apollo 11, 12 and 15 Basalts and Apollo IS Regollth) 
by reducing the nxm^er of trace elements Included. These trace element analyses 
are summarized In Table 6. Again, the open variance for SIO 2 was negative In 
the Apollo 15 regollth suite so that only three suites conform to the Chayes 
and Kruskal model. It Is apparent from Table 6 that all of the expected clo- 
sure correlations are positive and there is. In a practical sense, little 
difference between using one of these zero as the null value. From 

Equation (4) one can argue that the maximum positive value of p^^ can be 
approximated by*. 


ij 


o//C^Cj 


Equation (6) 


There appears to be a positive correlation between observed variance and hy- 
pothetical open variance for those arrays that conform to the Chayes and 
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Kruskal t^del and for most trace elen^nts 0 ^“ will be much greater than the 
Individual open variances so that the expression for maximum (Equation 4) 
can be used to examine the conditions chat must be satisfied to Induce posi- 
tive correlation. 

Results for a simple analysis of the requirements for positive correlation 
2 

(using the values for the Apollo 11 and 12 basalts) are given in Table 7. 

It is apparent from data presented in Table 7 that large positive closure 
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correlaclons are to be expected only between variables with a small product 

2 

of their coefficients of variation. As increases the product of the 
coefficients of variation required to produce a given maximum closure corre- 
lation Increases. Thus, unless the trace elements have very small coefficients 
of variation one need not expect to find intermediate in size positive closure 
correlations and the Chayes and Kruskal (1966) admonition concerning positive 
closure effects should be modified accordingly. As evidenced by comparing the 
data in Tables 2 and 1, relatively large positive null f irrelatlons are en- 
countered in the basalt and regollth suites but in all such situations they 
involve major elements with very small coefficients of varic.-tion. From 
Table 6 it Is apparent that the majority of the trace elements included in 
the three data arrays have coefficients of variation that are quite large. 

Shaw (1961) argued that when the coefficient of variation exceeds 0.20 to 0,25 
a lognormal distribution model was better than a normal distribution model 
for predictive and descriptive purposes 

A great deal of the chemical Information in the lunar data base consists 
of trace elen^nt analyses for which there are no accompanying major element 

analyses. In a strict sense such analyses do not conform to the Chaves and 

2 

Kruskal model as o^ can not be computed. In a practical sense, however, it 
seems appropriate to use the values of the coefficients of variation as a 
guide if one knows that all of the open variances for the major elements are 
positive for the suite from which the trace elements were drawn. If the 
product of the coefficients of variation is less than OrOl (an admittedly 
arbitrary selection) a positive null correlation of Intermediate level 
should be expected. If the product exceeds O.Ol and the Individual corre- 
lations are less than 0.25 then it appears appropriate to use zero as a null 
value. When the Individual coefficients of variation exceed 0.25 the required 
assumption of bivariate normality is probably not satisfied and alternative 


approaches are required to assess the degree of association between the 
two variables. 

Alternative Forms of Correlation 

In the proceeding sections the Pearson product moment correlation 
coefficient has been used as the measure of the linear association between 
a pair of variables « Tests of independence of this correlation coefficient 
require the assumption of bivariate normality which, as noted above, may 
not be met by geochemical variables. Demlrmcn (1976) recently reviewed 
the application of Spearman's and Kendall's rank correlation coefficients 
to situations in which the degree of association between a pair of variables 
is being examined. Both of these correlation coefficients are computed 
from ranks rather than the actual values of the variables and tests of inde- 
pendence are nonparanetrlc. As described by Demlrmcn (1976) the Spearman 
correlation ia obtained by treating the ranks as though they were the actual 
scores and is a measure of the monotonlclty between the two sets of scores. 
High values of Spearman's correlation indicate that the association is 
monotonlc (an increase in the rank of one variable is accompanied by an 
increase In the rank of the other or vice versa) although the relationship 
need not be linear. Kendall's correlation is a measure of the deerec to 
which the rankings of the two variables agree. 

Preliminary experlr^nts have been undertaken to examine the effect of 
closure of uncorrelated parent variables on the values of these two alterna- 
tive approaches to correlation, A set of 100 samples with 10 variables e.ich 
was drawn from uncorrelated parent variables of specified mean and variance. 
The computed Kendall and Spearman rank correlations all lead to ac^^eptance 
of the null hypothesis; that Is, the true values of tl;e correlations are 
zero and hence the ranks are independent. Ulten tills data set was clos< fy 
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transforming each row into percentage form and the Kendall and Spearman 
correlations computed for the closed set the null hypotheses are rejected 
for the association between variables with large variance. Thus, the 
preliminary work confirms intuition in that closure seriously modifies the 
rank correlations as well as the Pearson product moment correlation coefficient. 
There is apparently no procedure available for taking the effects of closure 
on the rank correlations into account. This is obviously an area that should 
be Investigated as the nonparametric aspects of hypothesis testing are cer- 
tainly appealing. 

Summary and Conclusions 

Many, but certainly not all, negative correlations of intermediate size 
are solely the result of closure. Similarly, it is possible that positive 
correlations of intermediate size arise as the result of the closure process - 
especially those between variables with very small coefficients of variation. 
Unless one makes use of the Chayes and Kruskal approximation method for 
assessing the effects of closure, however, recognition of the extent of the 
effects of closure and rational hypothesis testing are impossible. The 
automatic selection of zero as a null value is definitely to be avoided. 
Unfortunately, the practice of an Intuitive assessment of association of 
variation diagrams is widespread and almost always involves a picttire of ran- 
domness for which the Pearson product moment correlation coefficient is zero. 
Thus, the intuitive assessment amounts to a formal hypothesis test for inde- 
pendence using zero as the null but lacking in any quantitative input. 

From Table 2 and the accompanying discussion it is apparent that the 
effects of closure in the suites of Apollo 11, 12 and 15 basalt and regolith 
data suites is not negligible. Unfortunately, the set of Apollo 15 regolith 


samples does not conform to the Chayes and Kruskal model with the result 
that correlation coefficient based comparisons with other data sets is 
Impossible, It appears that there Is a very small probability that clo“ 
sure will Induce intermediate In size positive correlations between 
variables present In trace amounts In these data sets. Analyses of the 
trace element portion of a sample technically do not conform to the Chayes 
and Kruskal model but it has been suggested that a careful examination of 
the values of the coefficients of variation may allow the degree of associa- 
tion to be evaluated. Continued studies of the effects of closure on the 
Spearman and Kendall rank correlation are planned as the nonparametric 
nature of the tests for Independence is appealing for geochemical problems. 

The average compositions presented in Tables 1 and 6 may prove useful 
other investigators but a certain amount of prudence must be used as Is 
always the case when dealing with averages. 
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0.43 

0.079 

0.08 


0.47 

0.107 

0.04 

a 

tvj 

o 

0.06 

1.51 

0.04 


0.29 

0.346 

0.14 


0.08 

1.31 

0.04 


0.14 

0.112 

0.02 


0.20 

0.590 

0.20 

X 

to 

o 

0.09 

1.07 

0.04 


0.36 

0.307 

0.15 


O'. 08 
0.771 
0.01 


0.18 

0.521 

0.64 


0.13 

0.384 

0.04 

to 

o 

Cn 

o 

o • o 
• u • 
o o o 

O N> O' 
N3 


0.08 

0.391 

0.13 


0.05 

0.751 

0.00- 


0.11 

0.225 

0.04 


60*0 

ZZZ‘0 

61*0 

CA 

0.52 

0.263 

0.08 


0.35 

0.200 

0.06 


0.51 

0.314 

0.09 


0.31 

0.367 

0.87 


0.31 

0.249 

0.03 

n 

to 

O 

Co 


7.670 


29.720 


1.518 


6.719 


M 

CA 

ro 


15 . 


CA 




CA 

H 

M 

CA 

H 

H 

n 

CA 

*»3 

CA 

cS 

M 

H 

> 

O 
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W 


00 

> 

CA 

> 

CA 

\ 

00 


?3 
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o 

f 

H 


pa 


TABLE 2 


16 


COMPARISON OF OBSERVED AND CLOSURE CORRELATIONS* 


Suite 

Variables 

Observed 

Correlations 

Closure 

Correlations 

Apollo 11 . 

Basalts 

Si 02 -Al „03 

Si02“TiO2 

-.0129 

-.4583b 


-.3991a 

-.3057 


Si02-Fe0 

-.3532 

+. 3529b 


Si02“Na20 

-.3851 

+.2231 


AlnO^-FeO 

Ti02-Fe0 

-.5418a 

-.3431 


+.2828 

-.2626b 


FeO-MgO 

+.3014 

-.1457b 

Apollo 11 

Regolith 

Si02-Ti02 

- .4435a 

-.3836 


Si02-Fe0’ 

-.4977a 

-.3662 


Si02“Mg0 

-.5801a 

-.3043 

Apollo 12 

Basalts 

S 102 -A 1203 

+.2453 

-.1753b 


Si02“Pe0 

-.7147a 

+ . 3 /90b 


Si02“Mn0 

-.2326 

+. 4404b 


Si02-Mg0 

-.6297a 

-.7470 


MnO-MgO 

+.1487 

-.3589b 

Apollo 12 

Regolith 

Si 02 “Al 203 

+.0535 

-.4127b 


Si02-Fe0 

- . 5494a 

-.4593 


Al 203 “MgO 

-.5375a 

-.3005 

Apollo 15 

Basalts 

Si02-Ti02 

-.2867 

+ . 1340b 


Si02~Fe0 

- .6542a 

-.5130 


Si02-Mn0 

-.3269 

•<-. 3390b 


Si02-Ca0 

+. 5472a 

+.3284 


FeO-MgO 

+.2198 

-.3212b 


Fe 0 -Cr 203 

+.3018 

-.1019b 


Al 203 -Ca 0 

+.3223 

- . 1444b 

Apollo 15 

Regolith 

Open variance 

for Si 02 is negative 



* a indicates an observed correlation that is significant (99%) tested against 0 
b indicates an observed correlation that is significant (99%) tested against 
the closure correlation 
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TABLE 3 . 

CORRELATION COEFFICIENT MATRIX FOR APOLLO 11 SAMPLES* 
SiO^ AI2O3 Ti02 FeO MnO MgO CaO Na20 K2O P2O5 Cr203 S 


SIO2 1.00 + 

AI2O3 1.00 

TiO- 

FeO 

MnO 

MgO 

CaO 

Na20 

K2O + 

P2O5 

CT203 

S 


1.00 + + 

1.0 + 

+ 1.0 
+ 1.0 


1.0 


+ + 


+ + 

1.00 + + 

1.0 + + 

1.0 + 

1.0 


+ 

+ 

1.0 


TABLE ^ . 

CORRELATION COEFFICIENT MATRIX FOR APOLLO 12 SAMPLES 



Si02 

AI2O3 

TIO2 

FeO 

MnO 

MgO 

CaO 

N820 

K2O 

P205 


Si02 

1.00 

+ 



•• 


+ 




_ 

AI2O3 

4- 

1.00 

+ 

- 


- 

+ 

+ 

+ 

+ 

- 

TIO2 



1.00 



- 

+ 




- 

FeO 


- 


1.0 


+ 

- 

- 

- 

- 

+ 

MnO 

- 


+ 

+ 

1.0 

+ 






MgO 

- 



+ 

+ 

1.0 

- 




+ 

CaO 


+ 


- 

- 


1.0 

+ 



- 

Na20 



- 

- 

- 

- 


1.0 



- 

K2O 

+ 

+ 

- 

- 




+ 

1.0 



P2O5 



- 

- 




+ 

+ 

1.0 




- 


+ 



+ 




1.0 


1.0 

TABLE 5 . 

CORRELATION COEFFICIENT MATRIX FOR APOLLO 15 SAMPLES 
Si02 AI2O3 Ti02 FeO MnO MgO CaO Na20 K2O P2O5 Cr203 S 


Si02 l.OO + 

AI2O3 1.00 

Ti02 

FeO 

MnO 

MgO 

CaO 

Na^O 

K2O 

P2O5 

Cr203 

S 


+ 

+ 


1.00 + 

1.0 + 
1.0 


+ 

1.0 


+ 

+ 


1.0 


1,0 + + 
1.0 + 
1.0 


1.0 


1.0 


* + indicates positive correlation significant at the 997. level against closure, 
- indicates negative correlation significant at the 99% level against closure 


Y STATISTICS FOR TRACE ELEMENTS (PPM) IN APOLLO 11, 12, AND 15 BASALTS 
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TABLE 7. 

REQUIREMENTS FOR POSITIVE CU)SURE CORRELATIONS 


P£j (maximum) = 

0^2 « 21.54x10‘^* Oj. 2 »= 5.57x10"'^* 


Pj^j(max) 



pj^j(max) 

CiCj - Cj 

.80 

2.69x10-3 

.052 

.80 

6.96x10-4 .026 

.60 

3.59x10-3 

.066 

.60 

9.27x10-^ .035 

.40 

5.39k10"3 

.073 

.40 

1.39x10-3 .037 

.20 

1.08x10-2 

.104 

.20 

2.78x10'3 .053 


* Open variances scaled to agree with using the mean proportion Instead of the obser\ 
mean. 
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APPENDIX I. 

PRELIMINARY R MODE ANALYSIS OF APOLLO 15 BASALTS TAKING THE EFFECTS OF 
CLOSURE AND RATIO FORMATION INTO ACCOUNT 

Paper presented to the Conference on : Origins of Mare Basalts and their 
Inplicatlons for Lunar Evolution, November 17-19, 1975 
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* PRELIMIKARY R mDE ANALYSIS OF APOLLO IS BASALTS TAKING THE EFFIKTS 
CLOSURE AND RATIO FfNtMATION INTO ACCOU^n^: J.C. Butlir. Depirtmtnt of 
Coology. Unlvtrdty of Houiton. Noottoti, Tcxm 770<^ 

F«F«d with the Inttrpretetlon of < coUoetion of chtfmical an«lytes, oiw 
eon took for relationships between the variables (R node) or between the 
aampics <Q node). Ihe choice of a suitable measure of similarity for 
either Q or R mode ha been ditcuaaed by many authors and in this study, the 
Pearson product moment correlation coefficient is taken as the measure of 
similarity for an R mode analysis of selected Apollo 15 basalt chemical 
analyses as it 'provides a dimensionless index of the Joint behavior of each 
pair of variables with equal weighting of all variables. Each pair of 
variables can be examiited using a graphical (scatter diagram) or statisti- 
cal approach. Ihe goal of the examination is to determine which pairs of 
variables exhibit a significant relationship. To stake this determination a 
statement of what constitutes an insignificant relationship is necessary. 

In constructing a model to account .for the petrogenesis of a selected group 
of samples, one needs to be able to focus attention on those relations that 
significantly depart from randomness. If relationships that fail of signi- 
ficance are Included or if significant relationships are omitted, the model 
will not accurately portray the available information. Intuition and in- 
troductory statistics suggest that xeto Is the appropriate measure of ran- 
domness and if an observed correlation differs significantly from sero, 
the relationship is Judged to be significant. The purpose of this study 
» was to examine the chemical analyses of a suite of Apollo IS lunar basalts 
using three data forms in idiich tero may not be the appropriate null value: 
(1) closed data, (2) normalised date, (3) ratios with cMnion terms. 

Data secs in which the sum of the variables swasured for each sample 
is ci.nstant are termed closed and chemical data, by their very nature, are 
closed (1). Deviations between 100.00% and the sum of the measured variablea 
result from a combination of errors of commission and omission. It can be 
shown that there are no restrictions on the row bims of the variance-covari- 
ance matrix generated from open data (1). Open data being defined as data 
in which the constant item row sum restriction does not held. For closed 
data, however, each row of the variance-covariance utrix sums to sero: that 
is, Var^ + COV + COV gc’"*'*' 

Thus, for closed data, the sm of the covariances batmen A aiu) all other 
variables is equal to the negdtlve of vsria.,v* «f A which is positive by 
definition. There must be at least one negati covariance in each row of 
the variance-covariance matrix for closed data; there is no such restric- 
tion on the variance-covariance matrix for open data. As the correlation 
coefficient is simply the covariance divided by the product of the standard 
deviations, a negative covariance manifests itself as a negative correlation 
coafficient. In e data aat in which there waa sero correlation between the 
variablea in the open form, cloaure (by fomliqi pcrcentagea for example) 


/ 

/ 

/ 


R MODE AtOLlO IS BASALtS 


BUTLER. J.C. 


would Ifuiuct •one degree of correletton and, In general, the larger the 
virlaneea oi the larger the Induced correlation. Hiua, the obaerved corre* 
latton between varlablea In a closed data aet contatna a contribution due to 
cloauro aa well aa a Mature of the joint behavior of the variables. Chayet 
^1} argued that the contribution due to cloture should be used as the null 
value instead of aero. Ihls strategy often leads to a conflict with an in* 
tuitlve asseasnetit of randomness as large negative induced correlations are 
quite couBon; especially between the variables that are used In the fre- 
quently invoked Marker diagrams. Several siratigies have been suggested 
^ich nay appear to reduce if not elialMte the closure effect; (1) norne- 
lisation to a constant composition and (2) fomatlon of ratios, 

Sucss and Urey (2) appear to have intr^uced the strategy of normalit- 
tng chemical compositions to the composition of the average carbonaceous 
ehondrite by dividii^ the amount of each variable by the amount of that 
variable in average carbonaceous ehondrite (Cl). Ihus, comparlsion with 
average Cl is facilitated. Deamination of a normalised data matrix reveals 
that it is not doted in the sense that the row sums arc not constant. 
However, because each variable In each sample la divided by the same value 
(the ammint of that variable in Cl), the correlation between normalised 
variables (such as, SIO2 /SIO2 Cl versus CaO/CaOCl) is Identical to that 

between the raw or non-normallsed data (SIO2 versus CsO). 

Creation of ratios by dividing the value of one variable by the value 
of another variable In the same row does Indeed remove the effects of 
closure • a the row sums ere no longer constant. Depending upon the manner 
in which the ratios are formed, however, another form of induced correlation 
iBsy trite (s'tmmaflced In 1). For example. If the ratios have a common 
denominator (Ca0/Si02 versus AI2D3 /SIO2 ) there will be s correlation' 
induced between the ratios even though there may be sero correlation bet- 
iMcn the parents of the ratios, it is possible to compute and simulate 
the correlation betwtan ratios with common terms es functions of the co- 
efficients of variation of the parents. This expected correlation can be 
ueed as' the null value for determining whether the observed correlation 
differs significantly from that expected for rendomnisa. 

Chemical variations of Apollo 13 lunar basalts have been the subjects 
of many previous Investigations; Including, the recognition of olivine 
end quarts nromatlva subgroups. The apollo 13 basalts appear to be separable 
(on basis of major element veriatlons) from other lunar basalts (3) and 
this study is a precursor to e more detailed statistical analysis of major 
and minor elenent chemical variations of all of the lunar basalts. A total 
< of 28 major eleaient analyses of Apollo IS basalts were chosen as a trial 
•et of data (3 and 4). This aet of analyses was closed by recalculating 
Si02» AI2O3, tiO^i FeO, HgO, CeO, NSjO, K2O and p2^3^° lODX. The summary 
•tattatlea, varlanca-cevariaiiee lUtriw ami tba correlation coefficient 
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nttrix wtre computed end, to compare with previoue work, the correlattoni 
(r|j) were teated ageinat tero (Table 1 .). To uae the dtayea and Krunkal 
teat for atgniflcant departures from randomness one computes the means 
and variances of X - the hypothetical open array. Hits array is character- 
lied by having aero covariances atu) hence, tero correlations. When X is 
closed to form Y * the hypothetical closed array - the means and variances 
of the variahtes in Y are equal to those In the array being tested. As 
there are eeto correlations in X the correlations in Y (Pij) arise solely 
as a result of closing the data; these are the null values for the Chayes 
and Kruskal test. The open variances of the Apollo 15 lunar basalt suite 
arc all positive; a prerequisite for application of the test ( 1 ). Values 
of Pij must be considered as approximate as the equation for their calcu- 
lation is derived using a first order delta approx imst ion, A computer pro- 
gram In Fortran IV was written that enables sinxilstion of the X and Y arrays. 
The observed correlations (r^j) were also tasted against the simulated (qij) 
end epproximated (Pij) closure correlations and the results arc given in 
Table 1 . Of partieular Interest It that fact that correlations between 
Si02 - 1^.0 and AljOj - HgO would be Judged significant ( 9 V ?1 ^en tested 
against either p|j or qy. in fact, the expected closure correlation bet- 
ween SlOi” MgO would be significant ( 99 Z) if tested against rero. Yet, 
this la the correlation induced solely by the process of closing iltc data. 
Perhaps of greatest Interest, however, are the correlations between SiOi . 
Al 2 fb Fe(V;-igO vdtlch fail of significance when tested against zero yet, 
as the result of relatively large negative closure correlations, are 
significant when tested against both pjj and qjj. Thus, a perrogenetlc 
model for the Apollo 15 lunar basalts should include the correlations 
between Stn2 - At20j, and FeO - H| 0 . Utat these correlatlonb arc signifi- 
cant can only be detected if the effects of closure are considered. Simi- 
larly, such a model need not Include the correlation between SiOj - figO 
and AI2O3- MgO. 

Tha 28 Apollo 15 basalts were normalized by dividing each variable by 
the amount of that variable in average Cl (recalculated to 100'4 on a 
volatile-free basis; Mason). These basalts contain greater amounts of 
SIO2 , Al20^ Tin^, HnO, and CIO than the average cl and lesser amounts of 
Fee, MgO, N’ajo, K2O. and P2O5. Correlations that are r.tgnifleant against 
their appropriate null values are those given in Table 1 . 

Characterization and comparison of lunar sample chemistry by using 
ratios has been quite common. To Iflustrate the problems inhetent in 
the Interpretation of ratio correlations, the ratios Fe0/Al20j, FteO/Al203 
and Cao/Al203 were computed, if the parents of ratios with a comr>on 
denominator were uncorrelated and homogeneoua in coefficient of variation 
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t correlation of 0. S would be Induced between ration. 11tln In the corre* 
lotion referred to an Fcar-on (5). in Beneral, an the cocffl> ’ 

dent of variation o' the drntmlnator Inereancn relative to thone of the 
tmmeratora of the ration, the induced correlation increanen. Ratloi in 
t^ich the mifflcrators are In common or In wltlch the numerator of one ratio 
ia the denominator of the other ratio are also subject to a s purious 
correlation; that Is, a correlation that results solely from the formation ' 
of the ratios, observed (r^j), approximated (P^j), a^^^ simulated 

eorrelatloni between tba three ratios with M,.0^ as a common denominator 

A 3 ' . 

are given in Table 2. Tl»e parents of the ratios (see Table 1) are far 
from being uncorrelated yet large samples of retios formed from oncorre- 
lated variables with means and variances equal of those of the observed 
AlsOs. FeO, CaO and IfeO might exhibit correlations as strong as those 
exhibited by the ratios FeO/Al203* Can/Al202 and Ca0/Al203 - ^’.gO/Al2n3. 
there la a definite need for more experimentation with ratio formation 
before it will be possible to unainbigously assess the significance of 
ratios formed with common parts. Testing ratios with coenwn parts requires 
that the parent variablea be measured individually so that their means and 
•tarwlard deviations are known. Situations in which the ratio itself it 
measured instead of the individual variables are not testable. Thus, if 
one measures the ratios of Al/St and Ca/Sl, it will not be possible to 
isolate the component of the observed correlation that la the result of 
working with ratios with a common denominator tthe null portion) from the 
true correlation between the ratios. 

Vise of multivariate statistical models auch as cluster analysis, 
principal components analysis and various factor analysis models arc being 
used with increasing frequency In toe analysis of chemical data. All of • 
these models require a measure of similarity as a starting point. If the - 
similarity matrix is the correlation coefficient matrix or the variance- , 
covariance matrix then some method of taking the effects of closure 
normalication or ratio formation must be derived. If these effects are 
ignored, the results will be difficult to interpret as the investigator 
will not be able to separate out the bias that definitely ia introduced 
^en one deals with ratio; or percentage data forms. 
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TABLE 1. <*!**rvfd (tjj), Approxlmattd (Py) «nd Simutcttd (q^j) Corr«la- 

tloni for th« Apollo IS Lunar Batalt Data 


• 

variables 

r ^ 
"IJ 

••ij* 

91J^ 

SiOj-AljO., 

.1351 

-.2^21 

-.2951 

Sl02-Ti02 

-.7156 

-.0839 

.0733 

S102-F«0 

Sl02->te0 

-.7699 

. -.4409 

-.4' 32 

1 T 

3^ 

-.4712 

-.4710 

SiC2-CaO 

,5281 

-.1150 

.••.i2l7 

Al2<*3*foO 

Al20ymo 

Al203*CaO 

-.610A 

-.1634 

-.1272 

-.4517 

-.1634 

-.1836 

.6586 

-.0552 

-.0350 

Ti02*Fe0 

.8841 

-.0519 

-.0688 

Tl02-Mn0 

.4768 

.0498 

.0160 

FeO-MnO 

,4935 

.0049 

.0066 

FeO-tlfiO 

.3608 . 

-.2727 

-.^70 

FeO-Cae 

-.5336 

-.0822 

-,C942 

Mn(VKa20 

-.4242 

.1014 

' . 1202 

fig {> can 

-.9219 

-.1980 

-.1954 

N 820 -K 20 

.6833 

.0260 

.0383 

KjO* P'jO^ 

.5512 

,0099 

.0624 


1 an underlined rjj indicates significance against zero at 997. 

2 an underlined p^j Indicates significance against p^j at 9')?.. 

S an underlined q^^j indicates significance against qy at 997, 

TABCE 2. Observed (rij), Approximated (p^j^.and Simulated (qjj) Correia^ 
tions for Selected Ratios of the Apollo IS Lunar Basclt Data 


variables 




'’ti' 

FeO/Al 203 

CaO/AljOj 

.8488 

.7443 

, 7487 

Ie 0 /Al 203 

CaO/Al 203 

Mg0/Al203 

Hg 0 /Al 203 

.8260 

.5001 

.5164 

.5296 

.5002 

■ 5163 


1 an underlined rij indicates significance against zero at 997, 

2 an underlined p|j indicates aigniCleance against P|j at .997. 

3 an underlined qtj indicates significance against qij at 997. 
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APPENDIX II. 

VARIATIONS IN CHEMICAL C(»1P0SITI0NS OF APOLLO 15 MARE BASALTS 

Abstract of Paper Submitted to the Seventh Lunar Science Conference 
March 15-19, 1976, Lunar Science Institute, Houston, Texas 
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VARIATIONS IN CHEMICAL COMPOSITIONS OF APOLLO 15 MARE BASALTS; 

J.C, Butler, Department of Geology, University of Houston, Houston, Texas, 
77006. 

Sufficient major element analyses of lunar mare basalts have been 
accomplished so that within and between mission comparisons can be made, 

Hodges and Papike (1) have presented evidence that mare basalt texture is 
related to major element chemistry and that textural variations within each 
group of mare basalts is a function of cooling rate. Rhodes et al. (2) 

• effectively argue that careful investigation of major element variations may 
allow recognition of samples that retain information about their source area. 

Many investigators of lunar and terrestrial basalts resort to binary or 
ternary scatter diagrams (a form of R mode analysis) as the means of identi- 
fying groups of related samples and recognizing differentiation trends within 
each group. Ti 02 , for example, has been found to effectively discriminate 
against samples classified as to mission and plots of TIO 2 versus SIO 2 or 
KgO have been used to classify mare basalts (2). 

It has been demonstrated, however, that failure to take the effects of 
closure into account may lead to accepting correlations as significant when 
In fact they are probably due entirely to closure (3). In a set of 69 
chemical analyses of Apollo 15 marc basalts (Table 1) SIO 2 and MgO account 
for some 587. of the total variance of the data set but their observed 
correlation of -0.50 can be accounted for by closure effects (3). In addition 
reliance on scatter diagrams may result in a classification scheme In which 
the variables selected account for only a small amount of the total informa- 
tion content of the matrix. In the set of 69 Apollo 15 mare basalts (Table 1) 
used in this study, TIO 2 and MgO (used by Rhodes et al. (2) to group lunar 
basalts) account for only some 307. of the total variance. Sole reliance on 
R mode techniques for classification purposes may fall to detect relations 
between the samples. In fact, many methods used to relate samples on the 
basis of their chemistry are in fact discr inlnation techniques rather than 
true classification techniques. 

There arc a variety of techniques (0 mode), however, that are designed 
to classify samples as a function of values of the measured variables. Q mode 
cluster analysis and non linear mapping techniques were used to classify 
members of Che set of 69 analyses of Apollo 15 mare basalts. Interpretation 
of the results of such analyses is fraught with difficulty as a result of 
distortions in the sununary plots and pre-existing bias of the interpreter. 
However, the results of application of these classification techniques are 
interesting. As would be predicted the quartz normative Apollo 15 (pigeonite) 
basalts and the olivine normative Apollo 15 basalts are easily separated. 

In addition to these two groups it appears that there are two groups of 
Apollo 15 olivine normative basalts. Group 1 olivine normative basalts are 
characterized by low Si02 (66.86%), intermediate Ti02 (2.217.), high FeO 
(22.537.), high MgO (11.56%) and an MgO/CaO ratio greater than 1.0. Group 2 
^Hvlne normative basalts are characterized by intermediate SIO 2 (65.63%). 

Ti02 (2.567.), high FcO (22.53%), low MgO (9.62%) and an MgO/CaO ratio 
of less than 1.0. Values given In parentheses are group averages. 
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Normative composlt Ions? also emphasize the differences between groups 1 and 2. 
Using norms published by Rhodes and Hubbard (A) and Chappell and Green (5) the 
group 2 olivine normative basalts have an average of 10.27, normative olivine 
whereas the group 1 basalts have an average of 187. normative olivine (10 
and 11 samples respectively). The largest normative olivine content for a 
group 2 basalt is 14.7% and the smallest normative olivine content for a 
group 2 basalt is 16.6%; thus, the two postulated groups show no overlap in 
normative olivine content. Trace element data have not been used in this 
study but published values show some differences between groups 1 and 2. For 
example, the average Y content for type 1 is 28 ppm as compared with an average 
of 36 ppm (based on only 5 and 4 analyses respectively (5)). The same 
three-fold subdivsion (quartz normative and two olivine normative) can be 
recognized using non linear mapping and the recognized groups are identified 
by sample number in Table 1. 

Discriminant function analysis can be used to assess the assignment 
efficiency of the postualted groupings (although this involves admittedly 
somewhat circular reasoning). Si 02 and Ti 02 reclaim the postulated groupings 
with an efficiency of 94%. Nava (6) stated that the olivine-quartz normative 
dichotomy disappeared unless one used, the normative mineralogy or the Si 02 
content as variables. Q mode analysis was repeated after eliminating Si02 
as a variable with the result that the same three-fold subdivision was 
evident with MgO and FeO yielding an assignment efficiency of some 90%. 

It is interesting that the observed correlation of 0.24 between MgO and FeO is 
not significant (99%) if tested against a null of zero but is significant (997.) 
if tested against its expected closure correlation of -0.32 (3). Previous 
investigators may have overlooked the utility of the FeO-MgO plot to portray 
differences between samples on the basis that the observed correlation appeared 
to be insignificant. 

Previous investigators have noted that the olivine and quartz normative 
basalts can not be genetically related by a near surface crystal fractionation 
mechanism involving olivine, pyroxene or spinel (4,5). Variations within 
quartz normative Apollo 15 mare basalts, however, can be related to near 
surface olivine fractionation and cooling rate variation (1). Rhodes and 
Hubbard argued (4) that chemical variation within the olivine basalts was the 
result of olivine fractionation whereas Chappell and Green (5) believed that 
trace and minor element variations were indicative of derivation from more 
than one cooling unit. Along this line of reasoning, Papike and Hodges (7) 
reported that there are slower and faster cooling units or groups of Apollo 15 
quartz and olivine normative basalts. It is possible that the two sub-groups 
of the olivine normative basalts recognized in this study (Table 1) could 
reflect cooling rate variations although an in-depth interpretation of the 
significance of the postulated groupings is uncertain at the present time. 

In part, this reflects a lack of field control for the lunar basalts in the 
sense that there remains some uncertalnlty as to how many cooling units were 
sampled at the Apollo 15 site (4). 

There is no doubt that two groups of olivine normative basalts can be 
recognized from the Apollo 15 site using variations in major element chemistry 
as the basis for discrimination. All splits of the same parent sample arc 
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phncctl in the same sub>group and the classification reflects the views of 
previous investigators as to the relations between some of the smaller samples. 
For example^ Chappell and Green (5) considered that samples 15622 and 13636 
came from the same parent and that 15658, 15668 and 15674 were splits from 
the same parent. The former are assigned to group 1 and the latter to group 2 
(Table 1). The postulated two-fold subdivision of the olivine normative 
basalts may reflect: (1) a hiatus due to inadequate sampling; (2) the presence 
of two cooling units of olivine normative basalt at the Apollo 15 site; or (3) 
sampling of different parts of a single olivine normative basalt cooling unit. 

TABLE I 

CLASSIFICATION OF THE APOLLO 15 MARE BASALTS 

Group I » Olivine Normative - low SIO 2 

15016(3) 15636(1) 

15379(1) 15643(1) 

15555(6) 15659(1) 

15622(1) 15672(1) 

Group 2 - Olivine Normative - intermediate Si02 

15119(1) 15607(1) 

15256(2) 15658(1) 

15545(2) 15668(2) 

15556(1) 15674(1) 

15557(2) 15676(1) 

Group 3 - Quartz Normative - high Si02 

15058(1) 15495(1) 

15065(4) 15499(2) 

15076(3) 15595(1) 

15117(1) 15597(2) 

15475(2) 15118(1) 

15486(1) 
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mapping analysis in particular) to these data in an attempt to discover natural 
groupings based on alt of the major oxide data rather than relying on one or more 
of the many bintu'y scatter diagrams advocated by previous investigators. It 
appears that three chemically distinct groups of Apollo 15 mare basalts exist: (I) 
olivine normative basalts with tow SiO;, intermediate Ti02, and an MgO/CaO ratio 
greater than 1.0; (2) olivine normative basalts with intermediate SiOj, high TiOi, 
high MgO, and an MgO/CaO ratio less than 1 .00; (3) quartz normative basalts with 
high SiOj. low TiOj and FeO. Stepwise discriminant function analysis reclaims 
these three sub>groups with an assignment efficiency of greater than 95%. 


Apollo 15 Mare Basalts 

Sample making up the set of Apollo 15 mare basalts are identified as to the 
source reference, chemical, and textural properties in Table 1. Sample numbers 
given in Table 1 will be used to refer to the specific samples in the following 
illustrations. There are a considerable number of published partial analysis of the 
Apollo 15 mare basalts but only samples for which SiOj, AhOj, FeO, MnO, MgO, 
CaO, NajO, KjO, TiOj, and PjOj were analyzed were included. The set of 49 
analyses identified in Table 1 was selected to provide a training set of data; seven 
chemical analyses of Apollo 15 mare basalts (LSPET, 1972) were intentionally 
reserved as a set of unknowns in the event that a classification of the 49 training 
analyses proved feasible. 

Differences between the sum of the major oxides and 100% are the results of 
errors of commission and/or omission. Each analysis was recalculated to 100%. 
Closure to 100% changes slightly (0.01) the correlation coefficients but causes the 
sum of each row of the variance-covariance matrix to be zero. Summary statistics 
for the set of 49 analyses are given in Table 2. Si02, CaO, and FeO have very small 
coefficients of variation (.035, .058, .066, respectively) and hence are the least 
variable of the ten oxides. Of the oxides present in amounts greater than 1.0%, 
MgO and TiOi are the most variable (coefficients of variation of 0.18 and 0.17, 
respectively). 

The analysis of correlations between percentages of the same whole is 
difficult. There is no guarantee that zero is a suitable null value or that the same 
null value can be used for all possible pairs of correlations. It can be readily 
demonstrated (Chayes, 1971) that closure of uncorrelated parent variables will 
induce correlation as the sole result of percentage formation. Thus, the correlation 
between percentages contains a component that is a measure of the linear 
relationship between the variables. The problem is to ascertain the magnitude of 
the effect due to closure. 

Chayes and Kruskal (1966) developed a model that allows prediction of the 
closure effect. A hypothetical open array X (with zero covariances and correla- 
tions) is closed to form a hypothetical closed array Y, the variables of which have 
the same means and variances as those in the observed array U. Correlations in Y 
arise solely as the result of closure and thus may be used as null values. 
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Table 1. Make-up at the Apollo 15 mare basalt analyses. 


Sample 

reference Sampk 
number' numbe- 

CbemksI type 

PeUographic 

Source reference 

1 

155565 

olivine normative 

highly vesicular basalt 

Rhodes and Hubimrd, 1973 

2 

15256,15 

olivine normative 

recrystaUized basalt breccia 

Rhodes and Hubbard. 1973 

3 

15256.22 

olivine twrmative 

recrystaUized basalt breccia 

Rhodes and Hubbard, 1973 

4 

156685 

olivine normative 

highly vesicular basalt 

Rhodes aiul HuMuud. 1973 

5 

15545.13 

olivine normative 

porphyritic olivine basalt 

Rhodes and Hubbard, 1973 

6 

15016,4 

oliviiw normative 

hi^y vesicular basalt 

Rhodes and Hubbard, 1973 

7 

15555,8 

olivine normative 

porphyritic olivine basalt 

Rhodes and Hubbard, 1973 

8 

15475,35 

quartz normative 

porphyritic clinopyrosene basalt 

Rhodes and Hubbard, 1973 

9 

15499,2 

quartz normative 

porphyritic clinopyroxene basalt 

Rhodes and Hubbard, 1973 

10 

15058.4 

quartz normative 

porphyritic cUnopyroxene basalt 

Rhodes and Hubbard, 1973 

II 

150765 

quartz normative 

ptupbyritic cUnopyroxene basalt 

Rhodes and Hubbard, 1973 

12 

1507651 

quartz normative 

porphyritic cUnopyroxene basalt 

Rhodes and Hiridwrd, 1973 

13 

15118,4 

quartz normative 

porphyritic cUnopyroxene basalt 

Rhodes and Hubbard, 1973 

14 

15076,24 

quartz normative 

pmphyritic cUnopyroxene basalt 

Cuttitta etui., 1973 

15 

15065,31 

quartz normative 

1 

1 

n 

1 

1 

1 

Cttttitteetal., 1973 

16 

15486. lOA 

quartz normative 

7 

Cuttitta etui., 1973 

17 

15495,23 

quartz normative 

porphyritic cUnopyroxene basalt 

CuttitU etui.. 1973 

18 

15065.8 

quartz normative 

porphyritic clinopyroxene basalt 

Cuttituefoi.. 1973 

19 

15117,8 

quartz normative 

porphyritic cUnopyroxene basalt 

Cuttitta efai.. 1973 

20 

155573 

olivine normative 

porphyritic oUvine basalt 

Cuttitta efal., 1973 

21 

156075 

oUvine normative 

porphyritic olivine bualt 

Cuttitta etui., 1973 

22 

15659.4 

olivine normative 

porphyritic olivine basalt 

Cuttitta etui.. 1973 

23 

15555.27 

olivine normative 

porphyritic oUvine basalt 

Cuttitta etui., 1973 

24 

156435 

olivine normative 

porphyritic oUvine basalt 

Cuttituetul., 1973 

25 

15672,4 

olivine normative 

hi^ly vesicular basalt 

Cuttitta etui., 1973 

26 

153795 

olivine normative 

porphyritic oUvine basalt 

CuttittB etui.. 1973 

r 

15016.37 

olivine normative 

highly vesicular basalt 

Cuttitta et ol„ 1973 

28 

15676.6 

olivine normative 

porphyritic oUvine basalt 

Cuttitta etui.. 1973 

29 

15475 

quartz normative 

porphyritic clinopyroxene basalt 

Chappell sikI Green. 1973 

30 

15499 

quartz normative 

porph. clino, basalt vitrophyre 

Chappell and Green, 1973 

31 

15595 

quartz normative 

porph. clino. basalt vitrophyre 

Chappell and Green, 1973 

32 

15597 

quartz normative 

porph. clino. basalt vitrophyre 

Chappell and Green. 1973 

33 

15016 

olivine normative 

highly vesicular basalt 

Chappell and Green. 1973 

34 

15119 

olivine normative 

porphyritic olivine basalt 

CbappeU and Green, 1973 

35 

15545 

olivine normative 

porphyritic olivine basalt 

CbappeU and Green, 1973 

36 

15555 

olivine normative 

porphyritic olivine basalt 

Chappell and Green. 1973 

37 

15555 

olivine normative 

porphyritic olivine basalt 

Chappell and Green, 1973 

38 

15622 

olivine normative 

porphyritic oUvine basalt 

Chappell and Green, 1973 

39 

15636 

olivine normative 

porphyritic oUvine basalt 

Chappell and Green, 1973 

40 

15658 

olivine normative 

porphyritic olivine basalt 

(2happcll and Green, 1973 

41 

15668 

olivine normative 

highly vesicular basalt 

Chappell and Green, 1973 

42 

15674 

olivine nwmative 

porphyritic oUvine basalt 

Chappell and Green, 1973 

43 

15545.35 

olivine normative 

porphyritic olivine basalt 

Maxwell et al„ 1972 

44 

15555.153 

olivine mmnative 

porphyritic olivine basalt 

Maxwell etfll., 1972 

45 

15557.40 

olivine normative 

porphyritic olivine basalt 

Maxwell et of.. 1972 

46 

15065.6 

quartz normative 

porphyritic cUnopyroxene basalt 

Nava. 1974 

47 

15065,42 

quartz normative 

porphyritic clinopyroxene basalt 

Nava, 1974 

48 

15555,18 

olivine normative 

porphyritic olivine basalt 

Nava, 1974 


15597,19 

quartz normative 

porph. clino. basalt vitrophyre 

Nava, 1974 


'Used in the illustrations in the paper r ^fer to specific samples. 
Taken primarily from LSPET (1972) and i'hodes and Hubbard (1973). 



Vwiations in chemical composition of Apollo 15 mare basalts 


1433 


In order to obtain the expected closure correlations the means and variances 
of the X array must be assigned as described by Chayes (1971). Expected closure 
correlations can be approximated (p«, Chayes, 1971, Eq. 4.10) or simulated {q^) 
using Monte Carlo techniques. If negative open variances exist in the X array 
individual correlations can not be tested. This amounts to a crude rejection of the 
null hypothesis for the entire data set (Chayes, 1971) in that one or more of the 
correlations significantly departs from randomness but it is not possible to specify 
which one(s). 

Observed correlations (r^) were tested against zero, the approximated closure 
correlations ip^) and the simulated closure correlations (q^). Those that were 
found to be significant when tested against one or more of these null values are 
given in Table 3 along with their appropriate null values. 

That a number of correlations are significant when tested against all three null 
values is not particularly surprising. Of importance are correlations such as those 
between SiOr-AbOa. FeO-MgO, and FeO-MnO (Table 3) that have relatively 
small positive values and are not significant when tested against zero. Yet, these 
coirelations are significant when tested against their appropriate null values which 
are of intermediate size and negative. All three of these correlations would most 
likely be judged to be insignificant applying an intuitive appraisal to the scatter 
diagrams and would not be used in constructing a petrogenetic model. The 
FeO-M^ correlation, however, is extremely important in considering apparent 
sub-groups of the Apollo 15 mare basalts. 

Also of interest are the correlations between SiOj-MgO and AhOi-TiOj which 
are of intermediate size, negative, and significant when tested against zero. 
However, all fail of significance when tested against their expected closu:e 
correlations which are fairly large and negative. Such correlations should not be 
used in the development of a petrogenetic model except, perhaps, to consider why 
they were not significant. 

Rose et ai (1974) made use of what they referred to as correlation family 
diapams in which elements that were positively correlated were linked together. 
Signs of the significant correlations (against zero and against the ;;xpected closure 
correlations) are given in Table 3, and a schematic representation of the 
correlation families for those significant against the expected closure correlations 
is given in Fig. la. None of the oxides are in common with the two groups and. in 
general, there is a significant negative correlation between oxides in the two 
poups. What could be called the mafic group is made up of two sub-groups. 
Sub-group I links TiOi-FeO-MnO and could be considered as an opaque or oxide 
poup. Sub-poup 2 links FeO-MnO-MgO and could be considered as an 
olivine/pyroxene association. In the felsic poup, SiOr-Al.Or-CaO represents an 
anorthite association whereas P.O,-NajO-KjO suggests a KREEP component. 
This relationship (Fig. la) suggests that it would be possible to reduce the number 
of measured variables from the ten oxides to four (opaque, olivine/pyroxene, 
anorthite. and KREEP) and still retain a great deal of the information content of 
the data matrix. Thi.s would allow characterization of each sample in terms of the 
four new variaoles it contains. Procedures designed to accomplish this (such as 
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Table 3. Observe (r«), an>roxiiiiated (p«) and simulated ( 4 «) 
cotrelations. 


Variable I 

Variable 2 



9h 

SK), 

Al,0, 

0.26S 

-0.306 

-0.324 

SK), 

TiO, 

-0.6SJ 

-0.010 

-0.002 

SiO, 

FeO 

-0.77J 

-0.291 

-0.298 

SiO, 

MnO 

-0.36S 

0.193 

0.184 

SiO, 

MgO 

-OJJ0 

-0.516 

-0.510 

SiO, 

CaO 

0.492 

0.133 

0.132 

SK), 

Na,0 

0.314 

0.143 

0.142 

SiO, 

K,0 

0.40B 

0.042 

0.042 

Al,0, 

TiO, 

-0.319 

-0.065 

-0,039 

Al,0, 

FeO 

-0.464 

-0.231 

-0.200 

Al,0, 

Mi« 

-0.S82 

-0.085 

-0.073 

Ai,0, 

MgO 

-0.669 

-0.187 

-0.206 

AI,0, 

CaO 

0.370 

-0.147 

-0.132 

Al,0, 

Na,0 

0.630 

-0.033 

-0.022 

TiO, 

FeO 

0.869 

-0.062 

-0.079 

TiO, 

MnO 

o.sot 

0.041 

0.005 

FeO 

MnO 

0.495 

-0.017 

-0.016 

FcO 

MgO 

0.240 

-0.329 

-0.323 

FeO 

CaO 

-0.540 

-0.091 

-0.i09 

MnO 

MgO 

0.251 

-0.191 

-0.190 

MnO 

Na,0 

-0.462 

0.055 

0.067 

MnO 

K,0 

-0.299 

0.045 

0.032 

MnO 

P^>. 

-0.297 

0.046 

0.066 

MgO 

CaO 

-0.576 

-0.283 

-0.282 

MgO 

Na,0 

-0.621 

-0.069 

-0.080 

MgO 

K,0 

-0.5m 

-0.047 

-0.035 

MgO 

p,o, 

-0.3^ 

-0.052 

-0.073 

Na,0 

K,0 

0.641 

0.010 

0.001 

K,0 

p,o, 

0.649 

0.018 

0.043 

Na,0 

p,o, 

0.346 

0.001 

0.003 


'if fy is italic it is significant at the 99% level against a null of zero, 
’if Pm or (}y is italic the observed correlation (r^ ) is significant 
against the approprialc null value at the 99% level. 



R-mode principal components, cluster analysis, and factor analysis) are compli- 
cated by the closure effect. If, for example, one used zero :ts the null value against 
which to test the observed correlations, the coneiation families would appear as 
in Fig. lb, in which there is only the mafic opaque link. The olivine/pyroxene link 
disappears unless one uses the more appropriate correlation as the null value. The 
anorthite link is broken in Fig. Ib and, in addition. SiO] links with both NajO and 
KjO. If one were to use the standard computer routines for ^-mode analysis in 
which the matrix of correlation coefficients is used as the measure of similarity, 
the new variables would be similar to those in Fig. lb and the effect of closure 
would have been ignored. It has been suggested (Butler. 1976) that the matrix of 
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cjyiect^ cloture ccHTetatwns cmiM be subtrected from dw matrix of observe 
cmrelations and the difference matrix used as Um measure of similarity. Unless 
this or a similar strategy is lulopted, the effects of closure will contimw to influence 
the interpretation of the chemical data matrices. The danger is obvious. Without 
considering closure effects, one may both receive misinformation and lose valuable 
information. Preliminary results su^st that the proper ff-mode will be useful in 
characterizing the Apollo IS mare basalts and comparing them with mare basalts 
from other missmns. 


Classification of the Apollo 15 Mare Basali^ 

It is obvious that there are at least two chemically recognizable groups of 
Apollo 15 nune basalts— the olivine and quartz normative sub-groups rect^iz^ 
by many previous investi^tors. The SK), content, as might be expected, provides 
an efficient discriminant. In fact, Nava (1974) stat^ that unless SiO] w tlM 
normative amount of quartz were used that it would not be possible to recc^ize 
any groupings within the Apollo 15 mare basalts unless one considered coarse- 
grained basalts as a third group. Nava's statement is backed up with reference to 
an M^)-Ti 02 scatter diapam which does indeed reveal the difficulty in recc^iz- 
ing discrete groups. However, from Table 2. the pair of variables TiOj and MgO 
account for only 32.8% of the total variance of the 49 mare basalts — using 
Chayes's (1964) suggestion that the sum of the variances provkies a measure of 
the information content of the data matrix. It would appear advisable to base a 
classification scheme on the total data matrix if one is interested In ascertaining 
whether natural groupings exist within the data array. 

Scanning tlw raw data sug^sted that further subdivisbns of the Apollo 15 
mare basalts m^t be possible. For example, more than 90% of tlw olivine 
normative basalts contain either greater than 10% MgO and less CaO or greater 
than 10% CaO ami less MgO, 

There are several models that are designed to search out natural groupings of 
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samples tesed on variations of the measured variables. A number of these models 
were tried with Um hypothesis that a successful model would minimally reveal the 
olivine*quartz normative dichotomy. 

Q-mode cluster analysis is perhaps tlM simplest (conceptually) that could 
enable recc^ition of sample groups. However, results are difficult to interpret 
because of distortions that result in similar samples being placed in different 
groups or clusters. Results are depicted using a demlrivam whfeh is a tree*like 
plot in which pairs of similar samples are linked tether and treated as single 
samples. Groupti^ continues until all samples have been link^. A dendrogram for 
the 49 Apollo 15 mare basalts is given in Fig. 2 using the simple distance function 
(l^vis, 1971) as Um measure of between sample similarity. Any pre*existing 
of tt» investi^tor will certainly infliwnce the interpretation of such a dendro* 
gram. It is conceivable that a splitter would recognize the presence of 49 separate 
samples whereas a lumper would aigue for the presence of one poup composed 
of 49 entries. A similarity of i .20 was selected as the cut-off for group recognition 
because at this level tte olivii^-quartz normative dichotomy was preserved. 
Samples that linked below 1.^ were consideied to form a discrete p-oup of 
sub-set of the entire sampled populatkrn. Wlwn the cosine theta measure of 
similarity was used or when scaled or transformed (0.0- 1.0) data were used with 
the simple distance function the olivine-quartz normative dichotomy was not 
retained satisfactorily. 

Two major subsets of the sampled population are recognized in Fig. 2. Group 
A contains the oliviiw normative basalts and group B tlM quartz normative 
basalts. Within the olivine normative cluster (Group A) two sub-groups can be 
recognized: (I) A-l (Imnceforth referred to as p-oup (])) and (2) A-2 (henceforth 
referred to as poup (2)). The quartz normative basalts will be referred to as poup 
(3). 

As is freqiMntly the case in Q-mode cluster analysis there are samples that 
stand out as being sufficiently different fiom the major clusters. Sample number 23 
(15555,57) clusters with the combined poups I and 2 which suggests that. 
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although it is properly classified with the other olivine normative basalts, it does 
not fall into either of tiM postulatof sub^oups. Two splits of 15065 (sample 
numbers 18 and 46 in Table 1) arc .separated from the major clusters. These 
samples are characterized by iU>iMrmaliy low AI^Oj (5.36 and 6.05%, respectively) 
and high FeO (23.75 and 23.77%, respectively). Cuttitta et a/. (1973) suggested that 
pibbro 1 5065 protMibly contained two distinct ciMmical ami petro|p-aphic domains: 
(I) one with a relatively high mafic to plagioctase ratb (samples 18 and 46 in Table 
I) and (2) one with a mixture of mafic and felsic components (samples 31 and 42 in 
Table I). Chappell and Green (1973) considered that samples 15622 and 15636 
came from the same parent rock and that 15658, 156^, and 15674 also were splits 
from the same parent sample. The former are assigned as group I and the latter as 
group 2. 

In onier to assess the effects of possible distortion in Fig. 2, a model referred 
to as nondirwar mapping analysis (NLMA) was applied to the Apollo 15 mare 
basalt data. Howrath (1973) conducted a rather extensive series of tests of NLMA 
as applied to ^oic^ic data and a complete description of the method is given by 
Howrath (1973) and Salmon ( 1972). In general terms, the data are conskiered to be 
plottmJ in A# -dimensional space (where A# is the number of measured variables 
and the reference axes are orthogonal). NLMA attempts to locate a plane that can 
be passed through the cluster of points in Af -space upon which the projection of 
sample points reveals any natural groupings of points within the cluster with a 
minimum of distortion. The pair of vectors correspotKiing to the variables with 
neatest variance are used to define the starting projection plane and adjustments 
are made in an iterative fashion until the relationship between the data points 
projected onto the plane and the points in A# -space is maximized. In other words, 
the starting plane is tilted with respect to the selected vectors until the best 
possible two-dimensional representation of the Af -dimensional collection of 
points is attained. Experience with NLMA (Howrath. 1973 and Sutler, unpub- 
lished data) reveals that the groupings recognizable using NLMA are quite similar 
to those revealed by O’Utode form of analysis. Output from NLMA consists of a 
plot of the sample points on the final projection plane. It often occurs that the 
clusters of points would require greater than two dimensbns to reveal inherent 
arrangements within the data and an estimate of the probable dimensionality is 
given as part of the output. Inteipretation of the results of NLMA. however, is as 
difficult as with the (^-mode dendr<^ram as the investigator must make a decision 
as to how much space must occur between adjacent points before they are 
consbered as belonpng to different groups. 

The training set of 49 Apollo 15 mare basalts was treated using NLMA and the 
pbt of sample projection points is given in Fig, 3. The coordinate directions with 
greatest variance are those corresponding to the SiO- and FeO reference axes and 
the probable dimensionality of tiM plot is four. Thus, some distortion is present in 
the distribution of points in Fig. 3. As with the interpretation of the groupings in 
Fig. 2 there is a considerable degree of sub^ctivity associated with interpretation 
of the distribution of points in Fig. 3. Hie quartz normative basalts are well 
separated from the olivine normative group(s). Although it is possible to rec^nize 
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groups 1 and 2 of the olivine normative basalts one could reasonably argue that 
further subdivision is also possible. For example, the group 1 olivine basalts 
(located in the upper left hand corner of Fig. 3) could be subdivided using a 
boundary with a Y (the vertical axis) of value of 22.5. As might be expected, 
samples 18 and 46 appear again to form a distinct group and samples 13 and 23 are 
isolated from the major concentrations of points. 

In order to assess the hypothesis (Nava, 1974) that the olivine-quartz norma- 
tive dichotomy disappeared unless one used the normative mineralogy or the SiO: 
content of the samples, the set of 49 Apollo IS mare basalts was analyzed using 
NLMA without S 1 O 2 as a variable. Results of NLMA using the nine variables are 
given in Fig. 4. The probable dimensionality of the cluster of points is two and the 
coordinate directions with maximum variance are those corresponding to FeO and 
MgO. The probable dimensionality of two suggests that essentially the same 
information could be contained on a binary scatter diagram using FeO and MgO as 
coordinate axes. To facilitate further discussion, the postulated groups are 
outlined in Fig. 4. 
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Fig. 4. Non-linear mapping projection for the set of Apollo 15 mare basalts in which Si02 
was eliminated. Samples within a circled area belong to the groups defined in the text. 


Discrimination of the Apollo 15 Mare Basalts 

If one accepts the existence of three major groups of Apollo 15 mare basalts 
the strength of the classification can be assessed using discriminant function 
analysis (DFA). In general, most applications of DFA to chemical data (Chayes, 
1964; Chayes and Velde, 1965; Butler and Scotford, 1973) have involved situations 
in which the classification being examined was based on other than chemical data 
(such as modal mineralogy, age, structure, stratigraphic position, etc.). In the 
present investigation, chemistry is being used to assess a classification based on 
chemical data. Although admittedly circular, it was decided that this was 
worthwhile as a result of the uncertainty involved in interpreting the relationships 
in Figs. 2-4. If a classification based on the somewhat arbitrary classification 
techniques employed can be reclaimed with reasonable efficiency then one can 
begin to question whether such groups have petrologic significance. 

Samples belonging to the three groups are identified in Table 4. Samples 18 and 
46 were assigned to the quartz normative group and sample 23 was assigned to 
group 2 along with all other splits of 15555. Biomedical computer program BMDO 
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Table 4. Classification of the Apollo IS mare basalts. 


Sample number NASA 

Sample number (This report) 

Group 1— olivine normative — low SiOj 

13016 

6,27,33 

15379 

26 

15555 

7,23,36.37.44,48 

15622 

38 

15636 

39 

15643 

24 

15659 

22 

15672 

25 

Group 2— olivine normative — intermediate SiOz 

15119 

34 

15256 

2,3 

15545 

35.5,43 

15556 

1 

15557 

20,45 

15607 

21 

15658 

40 

15668 

41,4 

15674 

42 

15676 

28 

Group 3 — quaitz normative — high SiO, 

15058 

10 

15065 

15,18,46,47 

15076 

11,12,14 

15117 

19 

15475 

8,29 

15486 

16 

15495 

17 

15499 

30,9 

15595 

31 

15597 

32,49 

15118 

13 


M (stepwise DFA) was employed (Dixon, 1970 describes the routine and the 
method in considerable detail). Discriminant functions were computed for the 
three groups of Apollo 15 mare basalts in a stepwise fashion. The variable which 
had the largest influence on discriminating between the three groups was added to 
the discriminant functions at step 1 and additional variables were added to the 
functions using the same criterion. This continued until either all variables were 
added to the functions or until the addition of another variable did not change the 
discrimination. The stepwise procedure is highly appealing as the sequence of 
importance of the variables in the classification can be determined in addition to 
assessing differences between the groups. In the following discussion of assign- 
ment, efficiency is defined as the percentage of the total number of samples 
correctly assigned. 
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S 1 O 2 is the single most important discriminant. Ail 17 of the group 3, 13 of the 
group 2 , and 12 of the group 1 samples were correctly assigned for an assignment 
efficiency of 89%. T 1 O 2 was added at step number 2 with an increase in efficiency 
to 94%. In other words, the pair SiOj-TiO* provides a better discrimination of the 
49 samples into their three groups than any other pair of variables. The difference 
between this result and the choice of Si 02 and FeO as the coordinates of the 
NLMA projection plane results from the fact that DFA makes use of the within 
groups correlation coefficient matrix which effectively equally weights all vari- 
ables. In NLMA, on the other hand, the variables are weighted according to their 
variance. AI 2 O 9 is the third variable added which gives an assignment efficiency of 
98%. Addition of the remaining seven variables does not change this result. 
Coefficients of the three discriminant functions are given in Table 5. The 
discriminant functions are evaluated by multiplying the measured value of a 
variable for a particular sam>>lo by its appropriate coefficient and the function that 
yields the largest algebraic s- lution identifies the group to which the sample is 
assigned. Sample 27 (clas' .fied as a group 2 basalt. Table 4) is consistently 
assigned to the group 1 mare basalts as a result of its low (44.43%) SiOz content. 

DFA was repeated omitting Si 02 to assess Nava's (1974) implication that Si 02 
was not an essential part of a classification attempt for some coarse-grained mare 
basalts. Three variables lead to an assignment efficiency of 98%: (1) FeO (76%); 
(2) MgO (94%); and (3) AI 2 O 3 (98%). Coefficients of the three discriminant 
functions are given in Table 5b, 

Thus, the three groups of Apollo 15 mare basalts can be satisfactorily assigned 
using no more than three of the ten major oxides and recognition of the three 
groups does not require using Si 02 as a variable. Seven analyses of Apollo 15 
mare basalts (LSPET, 1972) were assigned to one of the three groups using the 
discriminant functions in Tables 5a and 5b. Sample numbers, chemical nature, 


Table Sa. Coefficients of discriminant function equations 
using ten oxides. 


Variable 

Group I 

Group 2 

Group 3 

SiO, 

537,409 

526.326 

559.876 

TiO, 

90.228 

87.091 

91.983 

AhO, 

784.926 

761.137 

796.706 

Constant 

-l.3607..<>66 

-13023.107 

-14703.089 

Table 5b. 

Coefficients of discriminant function equations 


using nine oxides (no SiO,). 


Variable 

Group 1 

Group 2 

Group 3 

FeO 

343.539 

351.576 

326.821 

MgO 

282.244 

289.326 

265.493 

Al,0, 

205.147 

213.315 

193.520 

Constant 

-5700.992 

-6017.839 

-5082.127 
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Table 6. Classification of the Apollo 15 mare basalts (LSPET, 1972). 


Sample number 

Clmnicsl t>-pe PeUograHiic type 

Classification 

13016 

olivine normative highly vesicular basalt 

Group 1 — tow SiO, 

15058 

Quartz normative porph. Uino. baaaJt 

Group 3— high SiO, 

15076 

quartz normative porph. dino. basalt 

Group 3 — high SiO, 

15236 

olivine normative recry. basalt breccia 

Group 2 — intermediate SiO, 

15499 

quartz normative porph. clino. basalt 

Group 3 — high SiO, 

15555 

olivine normative porph. olivine basalt 

Group 1 — low SiO, 

15556 

olivine normative highly vesicular basalt 

Group 2 — intermediate SiO, 


petrographic type, and assignments mre given in Table 6. Splits of the same sample 
are asHffted to die same p’oup and both sets of discriminant functions lead to the 
same assignment. 


Properties of the Three Groups of Apollo 15 Mare Basalts 

Textural and modal mineralogical data for the 49 member training set of 
analyses and the set of seven “unknowns” do not indicate a strong relationship 
between petrographic type (LSPET, 1972) and the three chemically defined 
groups. Of the five olivine basalt vitrophyres, three are classified as group 1 
basalts, and two as group 2. Ihere is a better agreement with Hava's postulated 
third group of coarse-grained olivine basalts and group 1 basalts although Nava 
(1974) was working with a fewer number of samples. 

Summary statistics for the three groups of basalts are given in Table 7. Group 1 
basalts have the lowest SiOi, intermediate TtOj, MgO greater than 10%, and an 
average MgO/CaO ratio of 1.22. FeO accounts for the greatest percentage of the 
variance and, for variables present in amounts greater than 1%, Ti 02 is the most 
variable. The hypothetical open variance of Si02 is negative with the result that 
individual correlations can not be tested for significant departures from random- 
ness. Group 2 basalts have intermediate SiOj, high TiO,, and an average 
MgO/CaO ratio of 0.?4. MgO accounts for the greatest percentile of the total 
variance and is the most variable of the major oxides. Again, the hypothetical 
open variance of SiO; is negative. Group 3 basalts are the quartz normative 
basalts with high Si02, low TiO,, and an MgO/CaO ratio of 0.85. The hypothetical 
open variance of SiO: is again negative. 

The normative compositions of the basalts also emphasize the differences 
between groups 1, 2, and 3. Using published norms by Rhodes and Hubbard (1973) 
and Chappell and Green (1973), the group 2 basalts contain an average of 10.2% 
normative olivine (ten analyses) whereas the group 1 basalts contain an average of 
18%. The largest normative olivine content for a group 2 basalt is 14.7% (15545, 
Chappell and Green, 1973) and the smallest olivine content for a group 1 basalt is 
16.6% (15555, Chappell and Green, 1973). Thus, there is no overlap between 
groups 1 and 2 on the basis of normative olivine. Similarly, normative ilmenite is 


Table 7. Summary statistics for groups 1. 2. and 3 Apollo 15 mare basalts. 

Group I 

SiO, Al,0, TiO, FeO MnO MgO CaO Na.O 
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high for the type 2 and low for the type 1. Trace element data have not been used 
in this study, but published values show some differences between groups 1 and 2. 
For example, the average Y content for the type 1 is 28 ppm as compared with an 
average of 36 ppm for type 2 (based on S and 4 analyses respectively (Chappell 
and Green, 1973). 

In all three groups, SiOi has the greatest mean and the smallest coefficient of 
variation. For closed variables with means less than 50%, the open variance will 
be negative if the square of the coefficient of variation is less than the sum of the 
open variances (Chayes, 1971). In effect, assessment of the strengths of linear 
relationships in the three groups is made impossible because SiOj is behaving as a 
constant within each of the three groups. In published analyses of the effects of 
closure on sets of chemical analyses (Chayes, 1971; Butler, 1975), the occurrence 
of negative open variances for SiOi has not been observed. However, published 
sets of chemical analysis usually reflect the desire of the investigator to illustrate 
some differentiation scheme and as a result an attempt is made to include a wide 
range of analyses. This effect is well illustrated in Table 2 giving the summary 
statistics for the combined set of analyses. The range of all variables is increased 
and there are no negative open variances. 

Various transformations have been designed to eliminate negative open 
variances (Chayes, 1971). However, these transformations have involved recalcu- 
lation of NasO to normative albite and/or recalculation of K 2 O to normative 
orthoclase with proper adjustments made to SiOj and AUOj. Because the total 
alkali content of the Apollo 15 mare basalts is very small the previously used 
transformations do not eliminate the negative open variances for S 1 O 2 in the 
individual groups. Until a suitable transformation is found, R-mode analyses of 
the Apollo 15 mare basalts must be confined to the complete set of data. 


Petrooenesis of the Apollo 15 Mare Basalts 

A considerable amount of information and speculation concerning the pet- 
rogenesis of lunar basalts in general and Apollo 15 mare basalts in particular has 
already been published and only a brief synopsis will be presented here. The 
group 3 quartz normative Apollo 15 mare basalts were collected from a wide range 
of geographic positions, samples were collected from Elbow and Dune Craters, 
the Apennine Front (Station 2), the edge of the Hadley Rille (Station 9) and at the 
ALSEP site. At Dune Crater, the quartz normative vitrophyres appear to underlie 
the porphyritic quartz normative basalts (Rhodes and Hubbard, 1973). Primarily 
on the basis of regolith compositions, Rhodes and Hubbard (1973) suggested that 
the olivine normative basalts were younger than the quartz normative basalts and 
that the quartz normative basalts were present at the surface as the result of 
excavation by impact. 

Most investigators have stated that the quartz and olivine normative types of 
Apollo 15 mare basalts can not be genetically related through a near-surface 
crystal fractional mechanism involving olivine, pyroxene, or spinel (Rhodes and 
Hubbard, 1973; Chappell and Green. 1973). Within the quartz normative group the 
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observed chemical variation can be related to varying amounts of fractionation of 
olivine at the surface or near the surface. 

Rhodes and Hubbard (1973) noted that although the majority of the olivine 
normative basalts were collected from Station 9A at the edge of Hadley Rille, 
there were samples collected that extended the range of this group to the 
extremities of the Apollo 15 site. The olivine normative basalts may have formed 
by a higher degree of partial melting than was involved in the genesis of the Apollo 
12 basalts (Chappell and Green, 1973). Chemical variations within the olivine 
normative basalts are related to near-surface olivine fractionation by Rhodes and 
Hubbard (1973) although Chappell and Green (1973) state that trace- and 
minor-element variations within the group are indicative of derivation from more 
than one unit. 

As a result of previous investigations of the genesis of the Apollo 15 mare 
basalts, it is necessary to question the validity of the two groups of olivine 
normative basalts postulated in this paper. There can be no doubt that the two 
groups can be very efficiently separated on the basis of modal and normative 
mineralogy, major-element chemistry, and perhaps minor-element chemistry as 
well. Does this reflect a hiatus to inadequate sampling or do these two groups 
represent two different lunar units? Analysis of the trace-element distributions 
may provide a means of differentiating the two groups. Occurrence of negative 
open variances of SiOi in the three groups of Apollo 15 mare basalts unfortunately 
precludes comparisons based on correlations between major and minor elements. 


Summary and Conclusions 

(1) If one uses the expected closure correlations as null values a set of 
significant correlations is obtained that differs from the set obtained using zero as 
the null value. 

(2) It may be possible to describe the Apollo 15 mare basalts in terms of an 
opaque, an olivine/pyroxene, an anorthite, and a KREEP component if significant 
correlations are identified using the expected correlations as null values. 

(3) Using Q-mode cluster analysis and non-linear mapping it is possible to 
recognize three groups of Apollo 15 mare basalts. Groups 1 and 2 belong to the 
previously recognized olivine normative basalt cluster and group 3 to the quartz 
normative cluster. 

(4) DFA reclaims the three-fold subdivision with a high assignment efficiency 
and the discriminant functions can be used to assign other Apollo 15 mare basalts 
to one of the three postulated groups. Only three variables are required to yield 
the most efficient discriminant functions. 

(5) Group I Apollo 15 basalts have low SiO^. intermediate TiO^, high FeO, 
high MgO, and an MgO/CaO ratio greater than 1 .0. Group 2 Apollo 15 mare basalts 
have intermediate Si02, high TiO:, high FeO, low MgO, and an MgO/CaO ratio of 
less than 1.0. 

(6) The apparent subdivision of the olivine normative basalts into groups 1 and 
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2 could reflect inadequate sampling at the Apollo IS site. If, however, the groups 
are real, in the sense that the chemical hiatus does indeed exist, then it is possible 
that the two groups represent different lunar units (flows?) or variations within the 
same unit with the group 1 basalts beitu{ enriched in normative olivine and group 2 
enriched in the non-olivine mafic silicate phases. It is also possible that the two 
groups are not genetically related but owe their genesis to partial melting at 
different levels in the lunar interior. 

(7) The occurrence of negative open variances for SiO; in all three groups 
precludes testing individual correlations for significant differences or for creating 
correlation family diagrams that might aid in elucidating differences between the 
two groups. Rational transformations designed to eliminate the negative open 
variance of SiOj are currently being investigated. 
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APPENDIX IV. 

Al/St VARIATIONS IN APOLLO 11, 12, AND 15 MARE BASALTS AND REGOLITH SAMPLES 

Paper presented to the Conference on O^mparisons of Mercury and the Moon, 
Houston, Texas, 15-17 November, 1976, 

A ccepted for publication In the Conference Proceedings. 
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Al/Sl VARIATION IN APOIXO 11, 12, AND 15 MARK MASAI.TS AND KKC.OMTII 
SAMTI.KS: John C. lluLlt-r, Ca*oloj*y, UnivorslLy ol liiuiston, Houston, TX 77004 

Tlu* iliroct mfiisuromi'nt of Al/Sl roll os in tho uppornufst portion of tlio 
lunar rrgollth by tlu* Apollo 15 and 16 missions makes one oploirlstic concern- 
ing possibillCies lor future unmanned geochemical surveys of the n.o«>n and 
other planetary bodies. These data have proven i-xtrimly usefiil In examining 
largo scale lateral chemical heterogenietlos and In mapping local chemical 
variations. Previous investigators have noted that the between site chemical 
variations reflect compositional differences in the parent material which Is 
primarly the local bedrock. The ability of variations in mi-arured Al/Si to 
mark the liighl and r. -mare material clearly has been demonstrated and a nearly 
complete trace and major element analysis has het-n obt.iined for the lunar 
highlands regolith by making certain assumptions and having access to prouml- 
truth observations. Wltliouc access to ground-truth chemical analyses (which 
will be the situation for initi.'il geochemical surveys of other planetary 
bodies) the present system will not yield absolute values of AI and SI. 
Chemical analyses of Apollo ll, 12 aivl 15 basalt and regolith samph's have 
b»*en extracted from the lunar data base* in an nttemiit to disct>ver how much 
Inhirmation in contained in the Al, Si, ami Al/Si variation (Table 1.) 

TAHI.K 1. 


S amp 1 e 



M«*nn SI 

Mean Al 

C Si 

C Al 

rAl :S1 

rAl :Al/Si 

Apol lo 

11 

basalts 

18.98 

5.00 

.02fa 

. 145 

-.044 

0.480 

A 1 ^ 0 1 1 o 

11 

Fines 

19.61 

7.21 

.02(* 

.016 

- .024 

0.782 

Apol lo 

12 

basalts 

20.87 

5.12 

.054 

.204 

. 248 

0.918 

Apol lo 

12 

Fines 

?1 .61 

7.10 

.081 

.081 

.412 

0.961 

Apol lo 

15 

basalts 

21.50 

4 . 84 

.04 7 

. 187 

.371 

0.954 

Apol lo 

15 

Fines 

21 .81 

7.70 

.019 

.211 

.379 

0.964 

Apol lo 

12 

and 1 5 

21.44 

5.91 

.044 

.294 

. 150 

0.977 


At a given sit«- the regolith samples tend ti> have higher Al, issenll.illy tin- 
same Si and higher Al/Si ratliis th.m the mare ha?alts frt»m the same 
The correlation between a ratio and its niimei alor is a simple (nnclUu) ol tin 
coefficients of variation of the parent variables (Al and Si) if the A1;S1 
tMiirelation is /.ero. As Si has a ri’laiively small coefficient «*f variation 
the correlations between A I and Al/Sl ai e quite large except for the Apollo 
11 lines. Plots of Al/Sl versus Al xhihit a strong positive linear trenil. 
There Is a consid<Tabl e difference betwicn the U'l’mls for the Apollo 11 
basalt and regolith samples. However, there is no statistical differince 
between the lines of organic correlation lor tlie Apollo 12 and Apollo 15 
basalt and r»*gollth samples. In fact, the Apollo 12 ami 15 basalt ami 
regiilith samples can be coiublned into a f ingle ilata set (Table 1) tor which 
the equation of organic correlation is Al*22.HD(Al /Si ) - 0.397. 'lln-se pre- 
liminary analyses stjggest that it may be possible to estimate the Al content 
of the lunar surface materials from the measured Al/Si ratios with relatively 
small assoclatid i-rri*rs. 


